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We show how to compute analytically time and space dependent correlations in one 
dimensional quantum integrable systems with an impurity. Our approach is based on a 
description of these systems in terms of massless scattering of quasiparticles. Correlators 
follow then from matrix elements of local operators between multiparticle states, the "mass- 
less form factors" . Although an infinite sum of these form factors has to be considered in 
principle, we find that for current, spin, and energy operators, only a few (typically two 
or three) are necessary to obtain an accuracy of more than 1%, for arbitrary coupling 
strength, that is all the way from short to large distances. As examples we compute, at 
zero temperature, the frequency dependent conductance in a Luttinger liquid with impu- 
rity, the spectral function in the double well problem of dissipative quantum mechanics 
and part of the space dependent succeptibility in the Kondo model . 
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1. Introduction. 

One dimensional quantum impurity problems arise in diverse areas of solid-state 
physics. Of recent interest are the tunneling through a point contact in the fractional 
quantum Hall effect [Q , the Kondo problem describing electrons interacting with dilute 
impurities in a metal, and the dynamics of double or multiwell systems coupled to an 
ohmic bath in the context of dissipative quantum mechanics , J| . A common feature to 
all these models is that they can be reduced to a model described by massless excitations 
in the bulk interacting with an impurity at the boundary. The absence of a mass gap 
leads to a power law behaviour for the current correlators in both the ultra-violet and 
the infra-red regime . The cross-over between these two regimes is non-trivial because of 
the renormalisation group flow induced by the impurity. A standard approach to these 
systems would be to use perturbation theory but it fails to capture all the physics, and 
new methods are needed. 

Beside (largely numerical) renormalization group approaches, another possibility for 
progress is provided by the integrability of some of these systems. Albeit not all of them 
are integrable, surprisingly, many are, and exact results can be obtained. The oldest such 
results concern the Toulouse limit of the Kondo model, which is a special point in 
the parameter space of the anisotropic model that is equivalent to free fermions. Other 
examples include the thermodynamic properties of the Kondo problem |J which were 
computed using the Bethe ansatz, and more recently the tunneling through a point contact 
in the v = 1/3 quantum Hall effect 0, where thermodynamic and as well some static 
transport properties were computed by a combination of Bethe ansatz and Boltzmann 
equations. Apart from the latter solution, transport properties as well as space and time- 
dependent properties are difficult to obtain because they require knowledge of correlation 
functions. 

Impurity problems of physical interest have a very simple bulk hamiltonian, typically 
a free boson. All the difficulty lies in the impurity interaction. In the basis where the 
bulk problem is simple this interaction is difficult to hanlde: in the classical limit, plane 
waves are scattered into very complicated wave packets by the impurity. If the theory is 
integrable, there is however another basis onto which the impurity has a simple effect: in 
the classical limit, there are some special wave packets that scatter simply at the impurity 
||. The price to pay of course is that in this new basis the bulk hamiltonian looks more 
complicated. Typically, this will lead us to describe a simple free boson as a gas of massless 
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quasiparticles (solitons, antisolitons and breathers) with non trivial scattering! The net 
result however is that both bulk and impurity are now manageable, and physical properties 
can be computed. 

Let us now recall that integrable massive theories, eg the bulk sine-Gordon model, can 
be described as a gas of quasiparticles with factorized scattering. A convenient formalism 
is to introduce creation and annihilation operators for these quasiparticles, like one would 
normally do for say a theory of free Fermions. If we denote by Z*(9) ( Z e {9) ) the creation 
(anihilation) operator of a quasiparticle of type e, the bulk interaction is encoded in the 
following Fadeev-Zamolodchikov relations : 

Z^(6 1 )Z^(6 2 ) = S e 7 e £ , 2 (0i - 02)^(02)^(01) 

t l t 2 

^(01)^(02) = 4,1(01 - 2 )^(0 2 )^(0l) (1-1) 

Z ei (0!)Z;(0 2 ) = 5 e e ^(0i - 6 2 )Z:, 2 Z<(6 1 ) + 27^5(01 - 2 )- 

In these expressions, the 0^'s are rapidity variables. In the discussion so far the theory 
is massive, and there is a meaningful mass scale m related to the coupling constant of 
some bulk perturbation. The energy and the momentum of quasiparticles can be then 
parametrized in the form E e = m e cosh(0) and P e = m e sinh(0), with m e oc m. For the 



relations (Ll) to make sense, constraints have to be satisfied by the S matrix, in particular 
it has to be a solution of the Yang-Baxter equation. 

It is then convenient to think of a massless theory as the limit of a massive theory 
when the bulk mass scale m goes to zero; for instance one can think of a free boson 
theory as the limit of the bulk sine-Gordon model when the amplitude of the bulk cosine 
perturbation vanishes^. Consider thus a massive theory, and take the following (massless) 
limit: = 0o ± (3 and m — » such that M = me 9 ° /2 remains finite. Excitations split into 
left and right movers, and dispersion relations become : 

E e = P e = M £ e^ rightmovers 

(1.2) 

E e = -P e = M c e~ p leftmovers. 

Here M is an arbitrary energy scale without physical meaning. The previous relations ( |1 . 1| ) 
are then decorated by a supplemental L or R subscript. The S matrices for the interaction 



1 Such a description requires in particular that the space of fields of the massive field theory 
considered as a perturbed conformal field theory and its massless ultraviolet limit are identical, 
which is the case in the "superrenormalizable" theories we consider here M. 
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between movers of the same chirality are unchanged, since in the initial massive theory 
they depended only on the rapidity difference §, but the right-left scattering becomes a 
constant phase (we refer the reader to H for a more detained discussion). This trivial 
left-right scattering will lead in the form factors approach to a left-right factorization, and 
will simplify matters drastically. 

This massless basis was not explicitely used in the original exact approaches to the 
Kondo problem. In the corresponding works, bare hamiltonians were used instead, and 
explicitely diagonalized via the Bethe ansatz. The point of view we take here is quite 
different: we work directly in the renormalized theory which is assumed to be factor izable 
||10|| (this following of course from the integrability of the bare hamiltonian) . This makes 
a tremendous difference when dealing with matrix elements of fields. While their compu- 
tation from the bare theory is extremely arduous, they can be obtained rather easily in 
the context of the factorized renormalized theory by using an axiomatic approach. The 



massless basis has proven very convenient in recent works ||11|| , in particular for what con- 
cerns static transport properties, like the DC conductance and the DC noise 0, [12] in 
the problem of tunneling between Hall edge states. Here, following the idea described in 
the pioneering work of [J13|, we show how this basis also permits the computation of time 



(or space) dependent properties. 

The natural approach to get correlations is to use matrix elements in the basis ( |1 . 1|) , 
the so-called form-factors. In the massive theory, these matrix elements are computed 
using a set of "axioms" ||14|| , ||15|| , ||16|| , generalizing Watson's theorem [I7|. Form- factors i 
of an operator O in a bulk theory are defined as : 

0nk,..., e „ =< 0|O(0, 0) (6>i) . . . Zl (6> n ) |0 > (1.3) 

where |0 > is the ground state, and their determination results form the following axioms : 
f(8i, 9i, 9 n ) ei ei>ei+1 €n Sj ' J +1 (9i - 9 i+ i) 



-f(9i, 6>i_|_i, 6>j, 9 n ) eij 



2 The S matrix depends on the relativistic invariant (E\ + E2) 2 — (P1+P2) 2 in the massive case, 
E1/E2 in the massless case, both functions of the rapidity difference using our parametrization 

3 These are actually called "generalized form-factors" since no order of the rapidities is 
prescribed 
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which is a consequence of (|LT|), and : 

/(0i, 9 n + 2m) ei ,..., en = f(9 n , B X) 9 n + 27ri) en>6l> ... )6ft _ is (1.5) 

which is a generalisation of the two particle form factor monodromy equations ||15|1 . The 
S matrix in the first relation has anihilation poles and this will induce similar poles in the 
form factors. Thus when 9 n = Q n -\ + in (which correspond to the mass thresholds in the 
two particles case for example) we have : 

i res /(6>i, n ) ei) ... :£n = f(9 u O n - 2 )e' 1 ,...,e' n _ 2 C en ^ n _ i x 
/ii i i \ (1-6) 

M c e ™-i e i (Q Q \ c r ™- 4e ™-3 (Q Q \n T »-3 £ n-2//J n \\ 

II — D Tl£l [Pn-1 —Ul)-— T „„ 3e „_ 3 [V n -1 — t/ n _3j^ e „_ l£ „_ 2 (,f n -l ~ I • 

Also, if the theory contains bound states, there are further relations involving the corre- 
sponding residues in the S matrix. For example, if there is a bound-state at 9 n = 9 n -i+iu, 
corresponding to a diagram depicted in figure 1. 




Fig. 1: Bound state diagram. 



Then, the corresponding relation for the form factor is : 

i res f(9 U 0n) eil ...,e n = a 7 (- 1 ) (aeB+1)/a ^n-i,6„/(»l. -> - ^k,..,e n _ 2 , 7 , (1.7) 



and here a 7 correspond to the square root of the corresponding residue in the S matrix 
for this process. 

Solving this set of equations is the most convenient method to obtain the matrix 
elements (although other methods based on lattice regularizations might also be applicable 
|TH] )- This is still a difficult task, but results are available for the sine-Gordon |15| and 



the sinh-Gordon [16] models, which is all what we need in this paper. Form-factors for the 
massless theory will be defined simply by taking the foregoing limit of massive form-factors 
(trying to formulate axioms for a massless theory per se might lead to some ambiguities) . 
Our strategy will then be to compute correlators simply by decomposing on intermediate 
states and using the exact matrix elements. This might sound a priori hopeless since there 
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are actually an infinite sum of relevant terms for operators of physical interest; however, 
we shall see that in many cases, the problem simplifies drastically. 

The paper is organized as follows. In section 2 we introduce the models and the 
quantities of physical interest which we want to compute. In section 3, the technique of 
form-factors is introduced with the example of the sinh-Gordon model. Although it is not 
of immediate physical relevance, this model is very simple since it has a single quasi-particle 
excitation, and appears useful pedagogically. In section 4, we introduce similarly the form- 
factors for the sine- Gordon model and compute the frequency dependent conductance in 
the v = 1/3 quantum Hall effect. In section 5, further applications are discussed; in 
particular we compute the spectral function in dissipative quantum mechanics and the 
uniform part of the space dependent succeptibility in the anisotropic Kondo model related 
to the screening-cloud problem. 

2. The models. 

The impurity problems we shall discuss here can all be mapped onto a hamiltonian of 
the form@ : 



where B is a problem dependent boundary interaction and the fields are defined on the 
negative half line. The geometry of the problem is shown in figure 2. 




(2.1) 
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Fig. 2: Geometry of the problem. 



4 In all what follows we set e = h = 1. 
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The method we present is quite general and in this section we proceed to an enumer- 
ation of the models we study and the corresponding quantities we want to compute : 



2.1. Sink- Gordon model. 

In this model, the boundary interaction is of the form : 

B = A cosh -(/>(> = 0,i). (2.2) 

Although we have no physical application for this model, it is much simpler than the 
models we will explore subsequently. The spectrum of massless excitations consist of only 
one particle and the form factors are much simpler than those used later in the sine- 
Gordon model. The quantities we will compute are the current-current correlation and the 
equivalent of the conductance in the Hall model. 

Our motivation for its study is to develop, in a simple fashion, the techniques used later 
in the real physical problems. We will show in a more lengthy manner how to obtain the 
current correlation functions and how the boundary affects these correlations. All features 
found in the other models are found in this one as well so it is a very good exercise. 

2.2. u = l/3 Hall effect. 

To start, let us recall the relation between the physical sine-Gordon model with an 
impurity and the problem in the half space. Let us start with the hamiltonian : 

1 Z" 00 1 

H = - / dx [8nvll 2 + —(d x ^) 2 ] + \S(x) cos(<^ L - <p R ), (2.3) 

2 J_ 00 Sttu 

where the L and R components depend on x, t as (pL(x+t), (pn(x—t). The bulk part of this 
hamiltonian physically describes the low energy edge degrees of freedom in the quantum 
Hall effect at filling fraction v = 2s 1 + 1 . When this Hall sample is subjected to a point 
contact constriction, transfer of fractionally charged excitations is possible. For generic 
filling fraction, many types of quasi-particles contributes as relevant charge transfer. The 
v = 1/3 is peculiar in that only the Q = e/3 charged laughlin quasi-particle is relevant 0. 
This is described by the impurity term in (|2.3| ). As discussed in @ in order to map this 
to a boundary problem, it is convenient to proceed in two steps. First introduce : 

(j) e (x + t) = [<pl(x, t) + vr(-x, t)} 

f (2-4) 

(p°(x + t) = — [<p L (x, t) - <pr(-x, t)] 
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which are both left moving. It is clear that the interaction term does not affect the even 
field, which therefore remains free. As for the odd term, it can be mapped onto a boundary 
problem as follows. Define : 

(f)° L (x,t) = Vty°(x + t), x<0, 

(2.5) 

<j)° R (x, t) = V2<j>°(-x + t), x < 0. 
The odd hamiltonian then reads : 

H=\( [8ng(U°) 2 + ^-{d x 4>°) 2 ] + XS(x) cos \<f>°, (2.6) 

and in the following we will write (p = <p° and g instead of v. Thus, for this problem, 
B — A cos \(j){x = 0, t). 

The quantity of interest in this case is the AC conductance at vanishing temperature. 
A standard way of representing it is through the Kubo formula : 

-j pL poo 

G(%) = - 5 pa / dx dx' / dy e ™ < j(x,y)j(x',0) >, (2.7) 

8tvuj m L 2 J_ L J_ oc 

where u>m is a Matsubara frequency, y is imaginary time, y = it. One gets back to real 
physical frequencies by letting um = —iw. In ( |2.7|) , j is the physical current in the unfolded 
system, j = d t ((fL — <Pr). Without impurity, the AC conductance of the Luttinger liquid 
is frequency independent, G = g. When adding the impurity, it becomes G = | + AG. 
After some simple manipulations using the folding, one finds : 



AG(um) = / dxdx' / dyc^ 

8nu M L 2 J_ L J_ 00 

[< d z ^(x,y)d^(x',0) > + < d s( p(x : y)d z 4(x' : 0) >} 



[2.1 



where z = x + iy. 

2.3. Dissipative quantum mechanics and Kondo model. 

The second model we will be interested in is the anisotropic Kondo model : 

1 r° 1 

H=- dx[8ngU 2 + (<9 x </>) 2 ] + A (S+e l ^ 0) / 2 + S-e' 1 ^' 2 ). (2.9) 

2 J-oo ^9 



It can be related to dissipative quantum mechanics || where it describes the dynamics of 
a double well system subjected to Ohmic dissipation. The hamiltonian for such a system 
is given by : 



hA he s 
H = — tt°x + iy°z + 2^ 



vl i , ( c, 



I I I u a ^ a \ ^QL n 

2m„ 2 \ rn a uj^ 



+ ^m a ul I X a ^-CT, 



(2.10) 



In this expression, the Pauli matrices cr x , a z act on the two dimensional space of states. A 
is a tunnelling matrix element and e denotes a bias between the two states. Dissipation 
comes from the coupling of this two states system to an environment of oscillators described 
by the last part of the hamiltonian. The effect of these oscillators (with mass, frequencies 
and coupling m a , to a , C a ) can be entirely encoded in the environment spectral function : 

It was shown [| that this system can be mapped, in the so called ohmic dissipation case, 
J (to) = 2ngu, to an anisotropic Kondo model. The conduction electrons in the latter play 
the role of dissipation and the z value of the spin is associated with the two states. In 
this paper we will work at zero bias (the bias would be equivalent to a magnetic field in 
( p.9|) ). The hopping A is related the the strength of the impurity A the precise relation 
being given in . This system has numerous physical applications |19| . 



The quantity of physical interest in the two state system is the effect of the bath on 
the tunneling between the two (degenerate minima). It is conveniently encoded into the 
correlator : 

C(t) = ±<[S z (t),S z (0)}>. (2.12) 

It describes the probability to be in a state S z (t) given that the system was in state ^(0) 
at t = 0. We will show how this is related to a current correlation and can therefore be 
addressed by the form factor approach. 

Another problem of interest is the screening cloud problem in the Kondo model. The 
three dimensional Kondo model can be reduced to one dimension by restricting to the s 
wave : 

^(r) = -4— [e~ ikFr Mr) ~ e ik ^ R (r)] + ■■■ (2.13) 
with the • • • denoting higher harmonics. In this language the Hamiltonian is given by : 

H = v F j dr ^\(r) l —^ L (r) - ^ R (r) % —^ R (r) (2.14) 



and the boundary interaction is 



B = XiP[(0)^ L (0)-S imp . (2.15) 

The screening cloud problem can be addressed by computing the local succeptibility 
which in 3 d is given by : 

X (r) =< ^(r)^(r) J dt S z tot >, (2.16) 



with : 



SLt = S? mp + I d\ ^(r)^V(r). (2.17) 



Here Si mp correspond to the impurity spin and the second part is the electron spin (the 
spin indices of the fermions are contracted with the Pauli matrix). When reducing these 



quantities to one dimension, we find following [2C| a uniform and 2k f part written : 



*( r ) = g^v bCun + 2 X 2k F cos(2/c F r)] (2.18) 



with : 



Xun =< [4«y^M+4Wy^Ml J dtS z tot > 
X2 kF =< [4(0y^W+^(r)y^(r)] J dtS~t 
and the total spin now given by : 

1 r°° rr it 

Stat = Sf mp + — dr [ViW-^(r) + ^(r)-^(r)]. (2.20) 

After having established the quantities we want to compute, we can bosonise|2l|]. Two 
bosonic fields are necessary: one associated with charge and one with spin. The charge 
field decouples completely and only the spin charge has interaction at the boundary. The 
action for the spin field is the of the form ( |2.9|) . The S z term in the action has been 
eliminated by a unitary rotation of the hamiltonian (which is unity at g = 1, the isotropic 
model). In terms of this bosonic field, the uniform part of the succeptibility at g = 1 takes 
the simple form : 

Xun(r) = \< d r <j>{r) j dt S z tot > (2.21) 
9 



(2.19) 



'tot > 



with nowu : 

Stot = S! mp + ^[ dxd^x). (2.22) 

We will show how to compute this uniform part at zero temperature. 

The problems are posed and our task is then to compute these correlation functions. 
This is a long story, and to explain the method we start by considering as a toy model the 
boundary sinh-gordon theory. 



3. Formalism: The sinh-Gordon model. 

3.1. The bulk current- current correlators. 

In most of this paper, we shall work in Euclidian space with x,y coordinates. Imag- 
inary time is at first considered as running along x. We consider the sinh-Gordon model 
with action : 

1 f°° r i 

S= / dxdy (d x <p) 2 + (c\0) 2 + Acosh</> . (3.1) 

lGTig J_ 00 L ' J 

This is a massive theory which is integrable; the conformal weights of the perturbing 
operator are h = h = —g. The spectrum is very simple and consists of a single particle of 
masse m, and S matrix : 

, N tanhi (d-iz£) , 

S{6) = 1) y., 3.2 

V ! tanhi + K ! 

where : 

2^7 



B 



1 + 9 

and 9 is the usual rapidity, P = msinh#, E = mcosYvO. Recall the duality of the S matrix 
in B — > 2 — B, i.e. in g — > 1/g. 

Consider the massless limit where A — » 0. In that limit, the current correlators are 
trivial and are given by : 

< d z 4>{z, z)d z '(f)(z\z ) >= — ^ 

< dz<p(z,z)8 M 4(z',z') >= - 

[z — z'Y 



5 This expression has to be considered with some caution but will be enough for the screening 
cloud computations of section 5. 
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On the other hand, we can formally describe this limit still using a scattering theory. 
As described in the introduction, this is done by writting 9 = ±(6*0 + 0) and taking 
simultaneously 9$ — > oo and m — > with me d °/2 — > M, M finite. In that case, the 
spectrum separates into Right and Left movers with respectively E = P = Me@ and 
E = —P = Me' 3 . The scattering of R and L movers is still given by ( |3.2|) where 9 — > 0. 
The RL and LR scattering becomes a simple phase, e Tl7rB / 2 . This phase will turn out to 
cancel out at the end of all computations, but is confusing to keep along. We just set it 
equal to unity in the following, that is we consider all L and R quantities as commuting. 

In this new description of the massless theory, we will need form factors in order to 
compute ( p.3| ). The form factors of the massive sinh- Gordon theory are well known ||22|| . 



We will only use the form factors of the fundamental field in what follows. By taking 



the massless limit of the formulas in |22[ , it is easy to check that can alter only the right 
or left content of states; in other words, matrix elements of between states which have 
different content both in the left and right sectors vanish. 

Our conventions are conveniently summarized by giving the one particle form factor 
of the sinh-Gordon field : 

< O|0(a;,y)|0 > R = ^exp [Me p {x + iy)] 

(3.4) 

< 0\(p(x, ?/) |/3 > L = ^exp [Me^x - iy)] . 

For the field 0, form factors with an even number of particles vanish. This is because as 
well as the creation operators of the sinh-gordon particle are odd under the Z 2 symmetry 
— > —(f). In the following we will use the notation : 

/(01, • • ■ , fon+l) =< O|0|ft, • • • , /W >R,...,R, (3-5) 

with the normalization of asymptotic states r < f3\(3' >r= 2ir5(f3 — (3'). Here, / depends 
on g, but we do not indicate it explicitely for simplicity. We have the following properties : 

< 010101, . . . , /Wl >R,...,R= (< O|0|0i, . . . , /Wl >L,...,L)* 

(3.6) 

< O|0|0 l5 . . . , 2n+ i >R,...,R=< O|0|0 2n+ i, . . . , 01 > L ,...,L, 



These form factors are expressed as follows 22 . Introduce : 



/mw(0) = Afexp 




dx sinh(^) S inh[^^] sinh(f ; 



X( ITT 



X 



sinh 2 (x) 



sm 



0)' 



2tt 



(3.7) 
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where : 



Then, 



J\f = exp 



-4 



dx sinh(^) S mh[^^l] sinh(f ) 
27 



£ sinh (x 



ft o a \ ( 4sin^ \ (2n+l) D / \ TT fminiPi ~ flj 

f(0 1 ,...,0 2n+ l) =fi\ 2 \ a V 2n+l 'P 2n +l(xi,. ..,X2n+l)[[ 



Fminii'K ' i -B) / ~\~ Xj 

(3.J 

where we introduced x = and the cr's are the basic symmetric polynomials : 



a p 



ii <i2 <---<i p 



with the convention ctq = 1 and a p = if p is greater than the number of variables. 
The P2n+i' s are symmetric polynomials, which can be obtained by solving LSZ recursion 
relations. The first ones read : 



P 3 (xi,...,x 3 ) =1 

P 5 (xx, . . . , x 5 ) =cr 2 a 3 - cfa 5 

P 7 (xi, . . . , x 7 ) =a 2 a 3 a 4 a 5 - cJ(<7 4 o| + a^a^e + a\a z - c\a 2 a b ) 
- c 2 ((Jia 6 a 7 + aia 2 a 4 a 7 ) + a 3 a 5 a 6 ) + c X (? 2 o^. 



(3.9) 



with ci = 2cos7rS/2, c 2 = 1 — c\. Observe that except for the overall normalization /j,(g), 
these expressions are invariant in the duality transformation g — > -. 

In ( p.4| ), fj, is an overall normalization for the form- factors. It is usually chosen by 
reference to the IR limit. However, we will require that the result ( |3.3| ) be recovered, and 
this involves a more complex computation. Using form factors, this two point function 
expands, assuming Re z < Re z' , as : 

< oi<w(,, wms, -010 >=-£/ (2 ')L;;(fr + +1 i)! M2 ( e/?1 + • • • + e/w ) 2 

exp [M(z - z') (e A + . . . + e^ 1 ))] • • • , /W)| 2 - 

(3.10) 

Now, by relativistic invariance, all the form factors depend only on differences of rapidities. 
Setting M(z — z') = e^°, (where (3q will in general be complex), one can shift all the /3's 
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by fio to factor out, for any 2n + 1 particle contributions, a factor ^ z ^ z ,y ■ Hence, the form 
factor expansion gives the result ([O]) provided iV is chosen such that 



J2hn+i=2g, (3.11) 



n=0 



where 



dpi . . . d/3 2n +i ( & , &„ +1 \2 (e ft + ... + A+i; 



= < (2*)»~(2n + 1)! (e " + ■ ■ ■ + e "' + " ' l/(ft - • ' 



|2 



(3.12) 

In practice, this sum cannot be computed analytically, but it can be easily evaluated 
numerically. The convergence is extremely fast with n, and for most practical purposes, 
the consideration of up to five particles is enough to get correct results up to 10 -4 . Similar 



convergence properties were observed in 

It must be emphasized that this result is very peculiar to the current operator. For 
most other chiral operators, the correct (z—z') dependence involves a non trivial anomalous 
dimension, instead of the naive engineering dimension. Hence, this dependence is not 
obtained term by term, as observed here, but rather once the whole series is summed up. 
Truncating the series to any finite n does not, in such cases, give reliable results all the 
way from short to large distances. The current is therefore an extremely favorable case, as 
would be the stress tensor, and we are fortunate it has a lot of physical meaning. 

3.2. Current current correlators with a boundary 

Having fixed the form-factors normalization, let us now consider the theory with a 
boundary. The geometry of the problem was illustrated in figure 2, where the boundary 
stands at x = and runs parallel to the y = it axis. The action is now : 

i f° [°° r 2 2 l f°° i 

S = / dx I dy (d x 4>) + {d y (p) +Acosh0 +A / dy cosh ~(f)(x = 0, y). 

Img J.oo J-oe L J J_ 00 2 

(3.13) 

This model is also integrable for any choice of A, A. The boundary dimension of the 
perturbing operator is x = —g. We can in particular take the limit where A — > while 
A remains finite. It then describes a theory which is conformal invariant in the bulk but 
has a boundary interaction that breaks this invariance and induces a flow from Neumann 
boundary conditions at small A to Dirichlet boundary conditions at large A. As discussed 
in H , the boundary interaction is characterized by an energy scale, which one can represent 
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as Tb = Me@ B . Tb is related with the bare coupling in the action ( |3.13j ) by A oc T B +9 . In 
the following, since obviously changes of M (which is not a physical scale) can be absorbed 
in rapidity shifts, we set M = 1. The effect of the boundary is then expressed by the 
reflection matrix : 



R{0) = tanh 



(3 TV 

2~ l A 



(3.14) 



In the picture where imaginary time is along x, the effect of the boundary is repre- 



sented by a boundary state. Following [24J] we can represent it in terms of the boundary 
scattering matrix : 



\B >= exp 



Ztv 



10 > . 



(3.15) 



In this formula, Z* denote the Zamolodchikov Fateev creation operators, K is related with 
the reflection matrix by : 



One can expand the boundary state into the convenient form : 



(3.16) 



\B>=J2 K([3b -Pi)-.. K{(3 b - • • • Z* L (f3 n ) Z R (Pi) ■ ■ ■ Z*(J3 n ). 

n=0 J0</3 1 <...<l3 n 

(3.17) 

Observe now that by analyticity, the matrix elements of d z (p between the ground state and 
any state with at least one L moving particle are identically zero. More generally, the only 
non vanishing matrix elements of d z (p are those where bra and ket have the same L moving 
part. The same results apply by exchanging d z with d z and L with R moving particles. 
As a result one gets immediately two of the four current correlators : 



<0\d z( p(z,z)d z ,(P(z\z')\0> 
< 0\ds<t>(z,z)d !i «t>(z , ,z , )\0 > 



2<7 



(z-z') 2 
(z-z') 2 



(3.18) 



which are identical with the ones without a boundary. 

The two other correlators are more difficult to get. Let us consider for instance : 



<O|a^,^)^0(z',z / )|5> 



(3.19) 
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The first non trivial contribution comes from the two particle term in the expansion of the 
boundary state : 



f°° dd 

/ JiK(p B -(3) < O\d- z 4>(z,z)d z ,0(z\z')Zi([3)Z* R ([3)\O > 

J-oo ^ 

f r fK(p B -P)e^exp[e' 3 (z + z' 

J-oo 27r 



(3.20) 



More generally, because \B > is a superposition of states with equal numbers of left and 
right moving particles, and d z (f>, respectively d z (f> act only on R, respectively L, particles, 
the expansion of (|3.19| ) takes a very simple form : 



£ 

n=0 



d(3i... d(3 2n +i 

( 27r )2n+l(2 n + l)l 



K(f3 B -fa). ..K(f3 B - 2n+1 ) (e^ +... + e^ 1 ) 



(3.21) 



exp [(z + z') (e ft + . . . + e^)] |/(A, . . . , f3 2n+1 ) | 2 . 

This correlation function depends on the product e^ B (z + z'). It is scale invariant in the 
UV and IR fixed point. These correspond respectively to sending (3 B to =Foo, that is 
the coupling A in the action to or oo, in other words Neumann or Dirichlet boundary 
conditions. In the first case, K = 1, in the second, K = — 1. Comparing with ( |3. 10|) and 
( |3. 1 1|) we find, as expected, that : 

2g 



< 0\d z (j)(z,z)d z ,<j)(z',z')\B >- 



(3.22) 



(z + z') 2 ' 

for Neumann, respectively Dirichlet boundary conditions. Although trivial, this result 
shows that the form factor expansion is well behaved, and allows us to study the correlator 
all the way from the UV to the IR fixed point when there is a boundary perturbation. 
In figures 3 and 4 we show the one particle (which is independent of B) and three parti- 
cles contributions. We observe that indeed the convergence, by looking at the respective 
contributions, is very rapid. 



0.03 



0.02 



0.01 



-0.01 



-0.02 



-0.03 




Fig. 3: One particle contribution. 
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Fig. 4: Three particles contribution for B = 1,0.1. 

The only drawback of this expansion is that it is not suited for studying the correlation 
of two operators right at the boundary. Indeed in that case, Re(z+z') = 0, and the integrals 
in ( |3.21[ ) do not converge. To solve this problem, we can introduce a modular transformed 
picture. We now consider the imaginary time as running along the y axis. Now the 
boundary is not represented as a state; rather, the whole space of states is different, since 
now we have only a half space to deal with. The asymptotic states are not pure L or R 
moving, but are mixtures. For instance, one particle states are : 

\\P>=\P> r +R(J3)\P>l- (3-23) 

More generally, asymptotic states are obtained by adding to . . . , (3 n >n l ,„ t R all combi- 
nations with different choices of R particles transformed into L particles, via action of the 
boundary. Only the following two terms contribute : 

P n >= \P U ...,P n >R...,R+... + R(j3 1 )...R(j3 n )\p n ,...,P 1 > L ,...,L +-... 

(3.24) 

Although we used the same notation as previously, different things are meant by L,R. To 
make it clear, we now use the conventions : 

< 0\(/)(x,y)\P > R = fiexp[e p (-y + ix)] 

(3.25) 

< 0\(j)(x,y)\p > L = /iexp[e /3 (-y - ix)]. 

To keep the notations as uniform as possible, we introduce the new coordinates : 

w(z) = iz = —y + ix, (3.26) 
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so here R movers depend on w, L movers on w. The normalization iV is of course the 
same as before, and as before the LL and RR correlators do not depend on the boundary 
interaction. One finds : 

< 0\d w (f)(w,w)d w >(f)(w' , w)\0 > = --. — 29 

(w — w'Y 

(3.27) 

< 0\du,(j)(w 7 w)d w x/>(w', w')\0 > = —f-z -to, 

(w — w'Y 



where we used the fact that \R{f3)\ 2 = 1. When compared with ( |3.18 ), these correlators 



have an overall minus sign due to the dimension h = 1, h = (resp. h = 0,h = 1) of the 
operators. 

Let us now consider : 

<0\diB<f>(w,w)d wf <f)(w' t w')\0> . (3.28) 

To compute it, we insert a complete set of states which are of the form (|3.23| ). In the 
massless case however, since d w (f> is a R operator, dwxfi a L operator, the only terms that 
contribute are in fact the ones with either all L or all R moving particles, as written in 
( ggg ). Thus, (JT2|) expands simply as : 

- E / \Jnir ^ -M- - 0b) («* + . . . + e*-) 2 „ om 

n=o J 1 > {6.2y) 
exp [(w - w')(e^ +... + e^)} \f(j3 u . . . , (3 2n+1 )\ 2 . 



Observe the crucial minus sign when compared to ( |3.21| ). It occurs because in one geometry 



the correlator depends on z + z' , while in the other on w — w' . This now converges provided 
y > y', even if x = x' = ie the operators are sitting right on the boundary. Now, using 
the fact that from factors depend only on differences of rapidities, this expression can be 
mapped with (|3.21| ) if we formally set j3 = (3' + provided one has : 

K(j3)=R(i~-l3), (3.30) 

as claimed above. 

To summarize, we can write the left right current current correlator in two possible 
ways. By using the boundary state one finds : 

POO 

< d s< />(x, y)d z ,<P(x\ y') >= / dE Q(E) exp [E(x + x') - %E{y - y')\ , (3.31) 
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(recall that x,x' < 0). One obtains Q(E) simply by fixing the energy to a particular value 
in (|3.21| ). When this is done, the remaining integrations occur on a finite domain for 
each of the individual particle energies since X^i" 1 = and there is no problem of 
convergence anymore. One then gets : 



dPt... dp 2n E< 



(27r) 2n+1 (2n + 1)! E - eft - . . . - eP*« 

(3.32) 



K(0 B - A) . . . K((3 B - (3 2n )K [p B - In (E - - . . . - e^ 2 ")] 
|/ [ft... /5 2m In (E-e^-...-e^)]\\ 



with the constraint Yli=i — The denominator might suggest some possible diver- 
gences; it is important however to realize that it vanishes if and only if the particle with 
rapidty P 2n +i has vanishing energy, in which case the form factor vanishes too. We can 
now shift the integrands to write equivalently : 



oo „o 

g ^ =E £-/-°° (27r) 2 "+ 1 (2n + l)! 1 - eft - . . . - eft™ 

K{\n{T B /E) - p x ) . . . K{\n{T B /E) - (3 2n )K [ln(T B /E) - In (l - - ... - e^ 2 ")] 
\f[P x ...p 2n M (l- e ^-...-e^)]| 2 , 

(3.33) 

where the constraint Y^i=i < 1 is implied, we used the fact that form-factors depend 
only on rapidity differences, and T B = e@ B . By using the dual picture, one finds : 

/•oo 

< d z 4>{x, y)d z ,<P(x', y') >= / dET{E) exp [-iE(x + x') - E(y - y')] , (3.34) 

Jo 



dpt... dp 2n 1 



where 



oo „o 

F(E) = -EJ2 / 

n=0 J ~ 



dP\ . . . dp 2l 



'2TT) 2n + 1 (2n + 1)! 1 - e& - ... - eft™ 
R(P 1 - \n{T B /E)) . ..R{P 2n - \n(T B /E))R [in (l - - ... - eft") - hx(T B /E)] 
\f[P 1 ...p 2n M(l-e^-...-e^)}\\ 

(3.35) 

where in ( 3.33 ) and ( |3.35| ) the constraint Yld=i < 1 is implied. The two expressions are 



in correspondence by the simple analytic continuation : 

Q{E) = %T{iE). (3.36) 
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3. 3. The analog of the conductance 

Although such a quantity does not have much physical meaning, we can formally 
define a conductance in the sinh-Gordon case using the current current correlators. It is 
instructive to carry out this computation now. 

To use the Kubo formula of the first section, we adopt the first point of view where 
the boundary is taken into account through the introduction of the boundary state \B >. 
Write again : 

POO 

<d^(x,y)d z '(j)(x , ,0)>= dEg(E)exp[E(x + x')-iEy]. (3.37) 

Jo 

This is the only correlation contributing to AG for positive Matsubara frequencies, and 

AG(u M ) = G -p^. (3.38) 

Here we have used the fact that uimL « 1, ie the system, although large, is much smaller 
than the wavelength associated with the (modulus of the) AC frequency. 

To go to real frequencies, we can simply substitute u>m — ► — in the K matrices in 
the integrals ( ^.33| ): 

AG(u) = —lm.Q(-iu) = -l-ReF(uj). (3.39) 
4o> 4a; 

Recall that K((3) = — tanh(/3/2). So we can expand the product of K matrices as a 
series, using 

K MTbM -0\= (fV - l) £ (3-40) 

Computing term by term gives AG as a power series in (uj/Tb) 2 - This is an IR expansion, 
valid for strong barriers Tb > oj. 

To get an UV expansion, holding for weak barriers Tb < a;, we have to split each 
integration into two pieces, Jl^^ 7 ^ and jli n ^/T B )- ^ n ^ ne ^ rs ^ integral, we expand K 
as in (|3.40|) but in the second case we expand it as 

K(ln(T B /iuj) — (3) = (l-I^e^) V (-TeY e -n(3 (341) 

This gives AG as a power series in (Tb/uj) 2 . 

A nice feature of the sinh-Gordon case is that the problem is well defined both for 
coupling g and for its dual |. This is because the dimension of the perturbing operator 
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being negative, it is always relevant, and none of these two couplings is plagued by problems 
of irrelevant perturbation theory, like what happens in the sine-Gordon case. We can 
therefore consider exactly the same problem with a coupling 1/g. On the other hand, these 
two cases are related by the sinh-Gordon duality, under which form factors are identical, 
up to an overall scale (due to the choice of \x as related with the two point function of the 
field </>). We deduce therefore the identity 

G L £-) = fG (3.42) 



T B J \g T b 

where the right hand term is computed using exactly the same formula as ( |3.21| ) and ( |3.38| ) 
but with the formal replacement g — > -. Tb = e^ B is the same in both cases. Of course, 
since the correspondance between Tb and the coupling A in ( |3.13j ) depends on g, ( |3.42| ) 



maps the conductance for g and 1/g with different boundary couplings. In other words, if 
we introduce the function T B {g, A) we have that 

% W ,A) = 5 2 G^, W ,A'j, (3.43) 



where A' follows from 



T B (g,\)=T B (±,\'y (3.44) 



Related duality properties are expected in the sine-Gordon model, but unfortunately are 
much more difficult to establish. 



4. v = 1/3 Hall effect. 

In this section we follow the same line of thought for the sine-Gordon model. This 
is the massive deformation of the free boson which preserves integrability with either 
boundary interactions B used in the Hall problem and the anisotropic Kondo model. Thus 
the form factors of the sine- Gordon model in the massless limit will be the quantities we 
need. The solitons/anti-solitons and breathers quasi-excitations make the problem more 
complicated but the results presented before hold with the addition of a few indices (and 
much more complicated form factors). 
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4-1. Expressions for the sine- Gordon form factors. 

We now consider the same problem in the sine-Gordon case, with action : 



1 

16ng 



S=^—l dxdy (d x (p) 2 + (d y( p) 2 + Acoscp . (4.1) 



To compare with standard normalizations, one has (3 2 = Sirg. The form factors approach 
is formally the same, albeit more complicated because the particle content is much richer, 
and depends on g. For 1/2 < g < 1, only solitons/anti-solitons appear in the spectrum of 
the theory This is the so called repulsive case, with g = 1/2 the Toulouse limit. When 
< g < 1/2, the particle content is enriched by [1/g — 2] bound states, called breathers. 
In the following we will denote by the indices e = ± the solitons and anti-solitons, and 
e = 1, 2, [1/g — 2] the breathers. The solitons form factors in the massive case were 
written by Smirnov and we obtain the massless form factors by taking the appropriate 
limit of the massive ones. Only right and left moving form factors survive in this limit, as 
in the sinh-Gordon case. Moreover, the symmetry of the action dictates that only form 
factors with total topological charge zero are non-zero for the current operator. As an 
example, the soliton/anti-soliton form factor is given by : 

<0\d x (t>(z,z)\P 1 ,P2>l£ R = 



e'^M(2nd) 2 e^ + ^' 2 n C( f 1 ^ exp [M(e^+e^)z] , (4 ' 2) 

cosh^GSi-ft + fcr) L J 

with e + e' = and e = ± stands for soliton (resp. antisoliton) . From []T3| one has : 



R r°° G in 2 *(fl+frQ Gin h *r(l-2g):c 

m = c sinh £ exp / . r * gx . r 2{1 ~ 9 ^ x dx \ , (4.3) 
2 I J i sinh 2(i- g ) sin h 7rx cosn ~Y 



with the constant c given by : 



a (-i u 1/4 / \ r oo sinh — sinh ^—^^ 
9 J I 4 Jo x sinh ^p^y cosh^ ?f 



and <i by : 



(4.5) 

and fj, is a normalization constant to be determined as before. Observe that in the free 
case g = 1/2, the form factors (f4.2| ) reduce to trivial kinetic terms since C(/3) oc sinh(/3/2). 

21 



The other soliton-antisoliton form factors follow from the analysis of [15|. Their 
expression simplifies in the case g = j, t an integer. This is the physically relevant case 
for the v = h fractional quantum Hall effect. One then finds : 



< 0\d z ^{z,z)\p 1 ,...,p 2n >r. 
t-1 



,-,+...,+. 

Ft. ...R 



sinh 



p=i 



p+n Hp 



(3 P - in) 



^M(87T 2 d) n e^ + - + ^/ 2 Y[ C(A - Pj) 

i<j 

n 2n 

\ | [ [ sinh _1 (t -l)(p q -p p ) detH. 

p=l q=n+l 



(4.6) 



The matrix H is obtained as follows. First introduce the function 

t-2 



ip(a) = 2 t ~ 2 Y[ sinh - ( a 



7TJ 7T 
Z h Z — 

t-1 4 



(4.7) 



One then defines the matrix elements as : 

„0 fc=2r 



Hij = — da TT V(a - /9 fc ) exp [(n - 2j - l)a + (n - 2i)(t - 1] 
2Z7T J_ 2i7r A = A 



(4.1 



where i, j run over 1, . . . , n — 1. It is not difficult to convince oneself that this produce 
a symmetric polynomial of the right degree. Although cumbersome, it is an easy task to 
extract these determinants, as examples we find for g = 1/3 : 



detH = exp 



1 2n \ 

Eft 



oi (e 



n = 2, 



2; i 



(4.9) 



cfetfT = exp(-^/3,) cr 1 (e /3p )(T3(e^), n = 3, 



up to irrelevant phases and with the cr 9 's defined previously. Having these expression we 
can get all form factors using the axioms they were constructed upon [JTS]]. For example, 
the solitons form factors with different positions of the indices we use the symmetry 
property : 

f(Pi, fa, Pi+i, 



Pn)e 1 ,...,e i ,£ i+1 ,...,e n S £ , ■ (Pi Pi+l) 



(4.10) 



—fyftli •••) Pi+l j Pit •••) /3n)ei,...,e i+ i,ei,...,e„- 

Here again, we omit the distinction between left and right moving form factor, they are 
simply related by complex conjugation. At the points g = 1/t the soliton S matrix used 
in the last expression is reflect ionless and basically just permutes the rapidities up to a 
phase. When there are breathers, the soliton S matrix has poles corresponding to the 
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bound states at the points f3 = m — 171 f° r the m'th breather. In view of the last 
relation, this induces poles in the form factors. We obtain the breather form factors from 
these poles : 



(4.11) 



res/(/3i,...,/3 n _i,/3 n ) ei ,..., e „_ lie „ = a m (-l) 2e " 2+lm C e „_ 1>e „ 

f(R R i l — _ 

J{Pl, Pn—l "T 2 Z 2^ _ ^ Jeiv, e «-2,7Tii 

and a m is given by the residue at (3 = in — ^ 9 g ^ m : 

a m = (res^ : -: : (/3)) 1/2 . (4.12) 

Having these relations, we posess all ingredients to compute all form factors for g = 1/t. 
Then, using them for the computation of the current correlations is merely an extension of 
the previous results for sinh-Gordon with indices. The normalisation of the form factors, 
H, is chosen such that (|3.18| ) is reproduced. This is fixed by introducing a complete basis 
of states : 



n=0 e 



[2-K) n n\ 



and computing the correlations exactly like in the sinh-Gordon case. It is interesting to 
observe however, that in the sine-Gordon case, there is another - a priori independent - way 
to fix the nromalization \i. Indeed, the operator d(j) being related with the U(l) charge, 
we need that 

/oo 
d x cf>\(3 2 >%= 2*6(0! -fa), (4.14) 
-oo 



using the fact that a soliton for the bulk theory (|4.1|) obeys </>(oo) — </>(— oo) = 2n. Using 
( fHf ) we get the requirement that 

"=S5cF*0 = (i^- (4 ' 15) 

Remarkably, this involves only the two particle form factor while the requirement that 
( p,18|) is obeyed involves a sum over an infinity of form-factors. However, the two should 
be identical if the description is consistent, which we checked is the case. 

Knowing the normalization before performing the sum ( |3.11|) gives us an a priori 
estimate of how many form-factors will be necessary to compute the full correlator. Indeed 
for g = 1/3, the one breather and 2 solitons form factors normalised to 3.14 which is very 
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close to the exact it. Similarly for g = 1/4 we found from the contributions up to two 
solitons that \i = 2.05 to compare with 2.094 = 2ir/3. 

Moreover the considerations concerning the correlations involving the boundary state 
follow in this case with now the boundary state given by : 



oo „ 

\B>=J2 K a ^((3 B -^)...K a ^{(3 B -^n)^2 ai (A)^ 6l (A)-^?"(^n), 

n=0 Jo<PK...<f3 n 

(4.16) 

with an implicit sum on the indices implied in this expression. The matrix K ab is related 
to the boundary R matrix in the following way : 

K ab ((3) = RI (4.17) 

The b means that we take the conjugate of the indices ie. ± — > =F and m — > m. 

From the previous expressions, we can compute de current-current correlation function 
in the presence of a boundary for g = 1/t. The results we will get depends on the boundary 
interaction, in the next subsection we present some general results for all values of g when 
the boundary is of the form ( |2.6[ ) . This is of relevance to the conductance in the quantum 
Hall effect. 

4-2. General remarks and analytic form of the conductance. 

To proceed we need the reflection matrix of the boundary sine- Gordon theory. For 

generic values of the coupling g, the amplitude for the processes H > + and > — is 

R±{(3 — /3b) , and for the processes H > — and ► + it is R^((3 — (3b) with : 

(l-g)i3 



R%{(3) =ie i ^r E R{(3) 



(4.18) 



where the function R reads : 



2g "4_ 
/ (1-9) _ 

4 g 2-rrg J " 1 V" 1 *- ) g 2-wg 

r(! + (' + ^)^- fl 3F)r(l + (; + i)^-^ 



r (I + ; H=»2 - r (i + (i + i)(i=sl - i<i=sM 



In ( 4.18 ), our conventions are such that in the UV limit {(3b — » — oo) the scattering is totally 



off-diagonal so a soliton bounces back as an anti-soliton, in agreement with classical limit 
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results for Neumann boundary conditions. A useful integral representation of R is given 
by : 



2 cosh 



2g 4 4 




dy . 2(1 -g)(3y smh(^)y 



R W = H^Z T ex P * 7e sin Z . : I ^7TT^7 • (4-19) 



2y £?tt sin h 2y cosh 



Recall that the spectrum is made of one breather and the pair soliton antisoliton in the 
whole domain 1/3 < g < 1/2. More breathers appear for g < 1/3, moreover the reflection 
matrix of the 1- breather is the same as in the sinh-Gordon case. There are no breathers 
for g > 1/2. 

In all these regimes, the form factors are known. They are quite complicated for 
generic g, and expressions for the correlators are more involved because the S matrix 
is non diagonal. We can however extract some features of the UV and IR expansions 
easily, following the same logic as in the sinh-Gordon case. To do so, consider the soliton 
antisolitons reflection matrix. Evaluating the integral in ( |4.19|) by the residues method 
leads to a double expansion of the elements R\, in powers of exp(/3) and exp(^ — . 
This leads for the conductance to a double power series in (uj/Tb)~ 2+2 ^ 9 and (uj/Tb) 2 
in the IR, (Tb/uj) 2 ~ 29 and (Tb/uj) 2 in the UV. Breathers do not change this result. For 
instance for the 1-breather, since the reflection matrix is the same as in the sinh-Gordon 
case, and therefore g independent, the contributions expand as a series in (uj/Tb) 2 in 
the IR, (Tb/uj) 2 in the UV. This holds for any coupling g. Therefore, as first argued by 
Guinea et al. J25|, at low frequency, the conductance goes as uj 2 for g < 1/2, u~ 2+2 / 9 for 
g > 1/2. The uj 2 power would seem to indicate that there should be a T 2 term in the DC 
conductance, but this is not correct because only the modulus square of R\, contribute to 
the DC conductance, and these expand only as powers of exp(^ — 1)P . 

The presence of analytic terms in the IR is a straightforward consequence of the fact 
that IR perturbation theory involves an infinity of counter-terms, in particular polynomials 
in derivatives of |[25|| . More surprising maybe is the fact that we find analytical terms 



in the UV. This requires some discussion. The UV terms follow from the short distance 
behaviour of the correlation function of the current. For any operator O we could write 
formally, 

00 r l i 

< 0(x>, y')0(x, y) >= ^(A) 2n dl...dn< d(x' , y')d(x, y) cos -0(1) ... cos -c/>(n) >, 

(4.20) 
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where O is the A — > limit of the field O. From ( |4.20| ), one would naively expect that 
the two point function of the current expands as a power series in A 2 , which would lead 
to a power series in (u /Tb) 29 ~ 2 ■ This is incorrect however because, even if integrals are 
convergent at short distance for g < 1/2, they are always divergent at large distances. It 
is known that these IR divergences give precisely rise to non analyticity in the coupling 
constant A. One usually writes : 

< 0(x', y')0(x, y) >= ]T C' 00 {x' -x,y'- y)Oi(x, y), (4.21) 

i 

where Oi are a complete set of local operators in the theory and the C"s are structure func- 
tions. These, being local quantities, have analytic behaviour in A. However, < Oi(x,y) > 
being non local is in general non analytic - actually, on dimensional grounds, 

<O l (x,y)>=A A /( 1 -^ocT^, (4.22) 

where A = h + h is the (bulk) dimension of the field Oi. If we computed the conductance 
perturbatively using Matsubara formula, we would use (|4.20| ) with O the electrical current 
operator. The case Oi the identity operator gives rise to an analytical expression in A, but 
eg the case Oi = d z dz4> gives A 2 /*- 1-5 -* times an analytical expression in A (its mean value 
can be non zero because there is a boundary) . More generally, since the only operators Oi 
appearing in the case of the current are polynomials in derivatives of 0, all with integer 
dimensions, we expect that the two point function of the current will expand as a double 
series of the form A 2n A 2m /( 1-9 ' 1 , ie going back to Tg variable , that the conductance will 
expand as a double series of the form (Tg / uj) 2n ^ l ~ 9 \TB /uj) 2m , in agreement with the form 
factors result. 

4-3. The free case. 

In the case g = 1/2 one has simply : 



R%{(3) = P{(3) 



el 3 



eP + i (4.23) 

In that case, only the soliton-antisoliton form factor is non zero, /(/?i,/?2) — /xe^ 1 / 2 ^^ 2 / 2 , 
where the normalization is easily evaluated \x = 2tc and we have set M = 1. Hence, ^(u) 
is readily evaluated 

/oo poo 
/ d^dl3 2 6(e^ + eh - u)[Q{(3 x )Q{(3 2 ) - P(/3 1 )P(/3 2 )]e^e^, (4.24) 
-oo J — oo 
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so 



x(l — x) + 1 
( x + iTf) (u-x + i^) 



Hu)=ooj dx , ■ ^ (4.25) 



from which it follows that 



Thus we find 



AG(u) = 1 - ^Un-\u/T B ). (4.26) 



G(u) = ~ (l - ^ tan-^w/Ts)^ . (4.27) 



This is in agreement with the solution of [26 . 



44. G(u) atg = 1/3. 

The conductance for (7 = 1/3 has a direct application to the quantum Hall effect. 
Comparing with the free case, previously treated, we now have a breather in the spectrum 
and non zero form factors for all number of rapidities. Still the convergence is such that 
evaluating the first few form factors give results to a very good accuracy, independently of 
the regime, UV or IR, in which we make the computation. 

In this case, the first few non zero form factors are /1, /±,=p, f±,^,i, f 1,1,1, etc... Here 
the subscript "1" denotes the breather. The first step is to compute the normalisation in 
order to satisfy ( |3.11|) . When computing this normalisation, we find that the first two form 
factors account for the whole result to more than one percent accuracy. Then including 
the 1 breather-2 solitons form factor is sufficient to get the result to a very good accuracy 
(-^1,1,1 is negligeable) . Actually one observes that the speed of convergence of the form 
factor expansion varies geometrically with the number of solitons (counting the breathers 
as 2 solitons). 

In order to get the conductance we need the reflection matrices, they were given in 
previous expressions and reduce to a simpler form for this value of g : 



1 r(3/8-^)r(5/8+ 2 ^ 

2 cos h(/3 - f ) r(5/8 - g)r(3/8 + g) 



R(p) = — 2 l J v ' ^ (4.28) 

y ! 2cosh(/? ' rr 1 — ■ " ■ 

and the breather reflection matrix is 



fl}G9)=tanh(^-^). (4.29) 
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From the pole of the 2 solitons form factor, the one breather form factor is found using 
( flll|) and its contributions to the conductance is : 



, Tvd 2 



log( 



-H z '^TZe tanh( " '^ Te 
8 2 



(4.30) 



here \x = 3.14 is fixed by ( |3.18| ) and d = 0.1414.... The contribution from the two solitons 
form factors is computed similarly, we find : 



AG(w) 



(2) 



^ -He I dp— — ; 2 T * \ ,f — Tb - \CW - log(l - e^)] | 2 



cosh 2 log (1 -e/ 3 )) 



e^(l-e^)(^ + 



(4.31) 

where C(/^) i s the function defined in ( [4.2|) . We can similarly write the following contribu- 
tion, and we find that these last two expressions are sufficient for any reasonable purpose, 
they give the frequency dependent conductance to more than one percent accuracy. We 
give the full function G(u) in figure 5. 
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Fig. 5: Frequency dependent conductance at T=0 



Observe that in the UV and in the IR we obtain the u dependance discussed previously, 
even with the truncation to a few form-factors. The form-factors expansion is indeed very 
different from the perturbative expansion in powers of the coupling constant in the UV, or 
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in powers of the inverse coupling constant in the IR. Each form-factor contribution has by 
itself the same analytical structure as the whole sum; contributions with higher number of 
particles simply determine coefficients to a greater accuracy. 



5. Anisotropic Kondo model and dissipative quantum mechanics. 

In section 2 we explained how the anisotropic model was related to dissipative quantum 
mechanics. The bosonised form of the hamiltonian is : 

H=\f dx Lyll 2 + ^-(d x 4>) 2 ] + $ (s + e^/ 2 + S_e-^°V 2 ) , (5.1) 
2 J-oo L ° n 9 J 2 V / 

where S are Pauli matrices. As for the quantum Hall problem, we keep using as a basis 
the massless excitations of the sine-Gordon model; however the boundary interation is 
different. This will result in different reflection matrices. Another difference with the 
previous section, will be the quantities we compute, the first step will be to relate them 
to current correlation functions and then using our by now well known techniques to get 
results. 

5.1. Dissipative quantum mecanics. 

We first work in imaginary time. We consider therefore the anisotropic Kondo problem 
at temperature T. Let us consider the quantity X(y) =< [S z (y) — ^(O)] 2 >. On the one 
hand, using that S z = ±1, it reads 2[1 — C(y)], where C(y) is the usual spin correlation 

C(y) = \[< S z (y)S*(0) > + < S*(0)S*(y) >] . (5.2) 

On the other hand, we can write a perturbative expansion for X(y) by expanding evolution 
operators in powers of the coupling constant A. At every order, we get ordered monomials 
which product of a monomial in and vertex operators of charge ±1/2. We 

must then evaluate S z (y) — S z (0) for each such term, trace over the two possible spin 
states, and average over the quantum field. Since we deal with spin 1/2, terms S + and 
must alternate, and there must be an overall equal number of S + and S~ , and an equal 
number of 1/2 and —1/2 electric charges. 

Now, since each S + (y) comes with a e -1 ^)/ 2 and each S~(y) comes with a e*^ y )/ 2 , 
S z (y) = S z (0) if there is a vanishing electric charge inserted between and y, and S z (y) = 
—S z (0) if the charge inserted between and y is non zero (and then it has to be ±1/2). 
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Therefore, we can write the perturbation expansion of X(y) in such a way that the spin 
contributions all disappear: 



i 00 z,t ry ry ry r^l J- r 

X(y)=-J2 x2n E £/ <*Vi / *!&•■■/ dy P / dy^... 

n=0 alternaUnqei=± p=0 J ° J ^ J Vp~^ J y J V 



4( 6l _|_ , . . _|_ £p ) 2 / e -itiHyi)/ 2 _ _ _ e -ie2 n 0(y2™)/2 



<%2n 

2/2»-l 

TV 

(5.3) 

where Z is the partition function, the factor 4 occurs because of the normalization S z = ±1, 
for every configuration of e's, only one value of ^(O) gives a non vanishing contribution. 
Here, the label iV indicates correlation functions for the free boson evaluated with Neumann 
boundary conditions (the conditions as A — > 0). 

On the other hand, let us consider the correlator 

< d x <f>{x, y)<t>(0, y') > N = -8g X (5.4) 

x 2 + {y- y'Y 

which goes to —8gTvd(y — y') as x — > 0. We have then, by Wick's theorem, 

< e -^4>(yi)/2 _ _ . e -«2„0(te„)/2 9a? 0( a . j y ) >N = 



2n 



X 



itiHyi)/ 2 p ^2 n 4>{y2 n / 2 )\ 



and therefore 

-» ei ^(vi)/2 _ _ _ e -i6 2 „0( 2/2n )/2 . y^g^x, y ') • 



(5.5) 



N 



2n \ / 2n 

X \ \ -< X 



jV 

i£2™</>(y2n)/2 



\~[ x +(y-Vi) Jytl + (y' - Vl ) 2 J \ 1 a 

(5.6) 

where contractions between the dots : are discarded. In ( |5.6| ), contractions between the 
dots would lead to a term factored out as the product of the two point function of d x (j) 
and the 2n point function of vertex operators, both evaluated with N boundary conditions. 
Now, we are going to be interested in the x — > limit where, with N boundary conditions, 
d x (f) vanishes. As a result we can actually forget the subtraction in ( |5.6|) , and write simply 
obtain 

X(y, A) = -7-^2 lim F F dy'dy" < d x cj>{x, y')d x ^(x, y") > A , (5.7) 

where the label A designates the correlator evaluated at coupling A, N corresponding to 
A = . Hence, we can get C{y) from the current current correlator. The latter can then 
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be obtained using form factors along the above lines. The only difference is the boundary 
matrix. If we restrict to the repulsive regime where the bulk spectrum contains only a 
soliton and an antisoliton, one has : 

Rl = tanh - ^ , R^ = 0. (5.8) 

Here again our conventions are such that a soliton bounces back as an antisoliton, in 
agreement with the UV and the IR limit that have Neumann boundary conditions. In the 
attractive regime we need to add the breathers with : 



= tanh 'f Sifflziii , (5.9) 



2 4(l/ g -l). 

tanh(| + 



Writing : 

poo 

< d- z cj>{x, y')d z cf>(x, y") > x = / G(E, (3 B ) exp [2Ex - iE(y' - y")] , (5.10) 

Jo 

we have that : 

lim < d x (f)(x,y')d x <f)(x,y") >\ = 

x^O 
poo 

/ dE [G{E, p B ) - G(E, -oo)] exp[-iE(y' - y")] + c.c, 
Jo 



(5.11) 



where the < d z (j)d z (/) > part and its complex conjugate (which are A independent) have 
been evaluated by requiring that the correlator vanishes as A — > due to N boundary 
conditions. Hence, using the fact that Q is real, 

1 f°° dl? 

X{y) = J o ^ [G(E, (3 B ) - Q(E, -oo)] sin 2 (£y/2). (5.12) 

Therefore, if we write : 

poo 

C(y) - 1 = / A(u) M ) cos(u M y) duj M , (5.13) 
Jo 



where um is a Matsubara frequency, we have : 

1 1 

An observation is now in order. From the foregoing results we see that 



A(u> M ) = , „_^ 2 [G(^>m,Pb) - G(uj m , -oo)] . (5.14) 



< 



1 r°° dF 

S z (0)S z (y)>-l = j^j o —ig(E,[3 B )-g(E,-oo)]cosEy. (5.15) 
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On the other hand, consider the expression 



/o i-O 
dx' / dx" [< d x <f>(x', y)d x <P(x", 0) > A - < d x cf>(x', y)d x <f>(x", 0) >,v] . (5.16) 
-oo J — oo 



By using the same representation ( |5.10| ), this is 

/0 r0 rOO 

dx' / dx" / dE [g(E, (3 B ) - G(E, -oo)] exp [E(x' + x") - iEy] + cc, 
-oo J — oo JO 

which coincides with ( |5.15|) after performing the integrations. We conclude that 

< S z (0)S z (y) > -1 =< J x (0)J x (y) >-< J x (0)J x (y) > N , (5.17) 
where we defined 

1 "° 



J * = ^— d x <j)(x,y)dx. (5.18) 
2^7r J_ oc 

We find also by the same manipulations that 

< S z (0)S*(y) > -1 =< J y (0)J y (y) >-< J y (0)J y (y) > N , (5.19) 

where 



1 



o 



J y = ^ d y <t>(x,y)dx. (5.20) 

2^7T J_ 00 

We now continue to real frequencies to find the response function : 

X» = i J dte^([S%t),S z (0)}), (5.21) 

to find : 

X » = -±—^Ll m [g(-i u ,p B ) - g(-iu, -oo)}. (5.22) 
{2gTi) z uj z 

As a first example, let us consider the so called Toulouse limit or free fermion case. 
Then the only contribution comes from the soliton antisoliton form factors, which as dis- 
cussed above is /(/?i,/?2) = [le^^e^l 2 . Hence, 

1 f°° r°° e^e^ 2 1 a 
x "(uj)= R e / dfixdfii—a 3-^7-5 jt^^ it^ Pi +e p2 -cA (5.23) 

that is 

X» = A^i m ( [" dx 1 ——) , (5.24) 

7t w Vio + iT b ){uj - x + %T B ) J 
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and 



4_Tb 

TV 2 OJ 



Im- 



1 



4Tl 



w + 2iT B 
1 



In 



w + iT B 



iT, 



B 



TV 2 ^ 2 + 4T| 



In 



T 2 



+ ttt - tan 
Tb 



Tb 



(5.25) 



In general, observe that, since the reflection matrix for solitons and antisolitons ex- 
pands as a series in e^ 3 , x ^ , will, for any coupling, expand as a series of the form (uj/Tb) 271 
in the IR. In particular, this leads to a behaviour C{t) oc ■h, t » 1 for any g. In the UV, 
one has to split integrals in two pieces as explained above in ( |3.41|) . Since the soliton- 



soliton form factors expansion involves powers of exp(- — 1)(3, 



expands double 



series in {T B /u) 2 ~ 29 and {T B /u) 2 in the UV. Hence at short times, C(t) - 1 oc t 2 ~ 2g '. This 



is in agreement with the qualitative analysis of [[27 



Results for g ^ 1/2 are more involved because there are non zero form factors at all 
levels. Still, when working out the first few form factors we observe a very rapid convergence 
with the number of rapidities and again we can give precise results for different values of 
g. As an example, let us show the results for g = 1/3. 

The computation for g = 1/3 is very similar to the previous conductance computa- 
tions. The boundary matrices are much more simpler though. In this case we have : 



tanh(- - — ), 



0, 



and : 



R i _ tanh (f " f 



tanh(| + f ) 



(5.26) 



(5.27) 



Then, as was found for the conductance, we find that the first two contributions are 
sufficient for most purposes, they are given by : 



9fi 2 d 2 
8tvuj 



Tie 



tanh 



log( 



V2T, 



-) 



tanh ( + ^ 



(5.28) 



and : 



5x» (2) 



( 3\xd 



-Tie 



to 



d(3 



|C(l-log(l-e^))| 2 
cosh 2 (/3-log(l -eP)) 



[R±(/3 + log(w/T B ))i2i (log[(l - e?)uj/T B ]) - l] . 

(5.29) 
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Fig. 6: Spectral function for Tg = 0.1. 



w 



Again these two expressions are sufficient to get a very precise result. Similar com- 
putations give rise to the results in figure 6 where we plotted x"(u;)/u; for the values, 
= 3/5,1/2,1/3,1/4. 

When making these calculations we have to be careful about which terms are needed 
for a good convergence. Our observation is that keeping the form factors up to two ra- 
pidities give very good results precise to 1%. It is possible to go further and get a better 
precision if needed. The last statements are true for g G [0.6, 0.2] and we believe even 
further (we can get rough bounds on the higher contributions and have an idea of the 
precision). Still, for the moment the isotropic Kondo point is difficult to treat. 

A surprising result is that the emergence of quasi-particle peaks in S(u) is found at 
= 1/3 and not at g = 1/2 as was expected from other means of calculations. Physically 
this means that the behaviour of the two state system goes from coherent to incoherent 
behaviour at that value of g. This is supported by a recent RG numerical study |?8 . 



5.2. Shiba's Relation. 

Up untill now, we showed results for certain values of g more or less limited by our 
ability (or tenacity) to write the form factors corresponding to that value of the anisotropy, 
and make them converge. It is not impossible to find general relations though; for example 
the behaviours in the UV and the IR in different models were infered in all generality. 
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Here we present a generalisation of Shiba's relation [S9J which was proven for the An- 



derson model and generalised to Luttinger liquids by Sassetti and Weiss|30]]. The relation 
states that : 



lim 



with xo the static succeptibility. If we look at the quantity : 



(5.30) 



oo „ 
. n J — C 



d(3i . . . d(3 2r 



1 



2n) 2n + 1 (2n + 1)! 1 - eP* - ... - 
K a ^(ln(T B /E) - 1 ) . ..K a — b ^{\n(T B /E) - n _ t ) 
K a nbn [i n (r B /E) - In (1 - - ... - e^ 2 ")] |/ [A . . . fo n , In (l - 



n=0 ' 
ai6i , 



(5.31) 

insert it in the expression for x"(w) an d expand it around E ~ we find that the contri- 
butions from the K matrices all cancel and only a constant is left (we have to take into 
account the fact that the soliton/anti-solitons K matrices always appear in pair). Then 
comparing this with the UV normalisation we find that : 



)] 



lim 



7T 



(5.32) 



The total succeptibility is x = x' + ^x" an d the static succeptibility xo which is the zero 
frequency limit of x' can a l so be infered from the previous expressions for the spin-spin 
correlation. We just need to take the real part when continuing ( |5.14| ) to real frequencies, 
which leads to : 

1 

Xo = ■ (5.33) 



Tv 2 gTj 



B 



Then in order to make contact with the previous expression we need to renormalise the 
spins to 1/2 and put the correct normalisation!!. This amounts to multiplying the each of 
the previous expressions by tv/2 which leads to the correct result : 



lim*>> 

UJ— >0 UJ 



(5.34) 



5.3. Screening Cloud problem. 



Another problem that can be addressed using the form factor techniques is the screen- 
ing cloud problem for the anisotropic Kondo model. This has been a long standing problem 



6 We compute |[»S Z , S z ] but the succeptibility has ^- in front instead of 1/2 thus at h = 1 we 
have to renormalise by 2tt. 



35 



and many theoretical studies were devoted to it Ell] . Our one dimensional formulation here 



follows from [20 



Our aim here is to compute the uniform part of the succeptibility, defined in the 
introduction. We restrict to the case g = 1 in what follows. The electron spin density is 
given by : 

1 r° 

d x 4>. (5.35) 



2 

We can proceed to the evaluation of Xun, which is given by : 

A Xun(r) = -^Jdt< d r (/>(r)Stot(t) > ■ (5.36) 

In the following, we will always substract the free part (Tb = 0, (3b = — oo) and denote the 
corresponding quantity by the symbol A. To find correct results for the zero temperature 
succeptibilities, we always need to take the length of the system to infinity first before 
taking T — > 0. In this case, this means that we always do the time integrals last. 

Let us compute separately the static electron-electron and electron-impurity, Ax ee , 
Axei succeptibilities. The electron-electron contribution is the appropriate integral of the 
correlator < d x <fid x (f) > — < d x (j)d x 4> > n '■ 

A Xee (r,T = 0) = — / dy dx dE [Q(E, T B ) - G(E, 0)] e E{x+r ^ y + c.c, 

= j_ [g(E,T B )-g(E,o)] cEr 

8tt e^o E 
= 0, r ^ 0. 

(5.37) 

with y the imaginary time. For r = 0, the integral over E is divergent and this computation 
does not make sense. In fact, Ax ee (r, T = 0) has a delta like contribution at r = as we 
now show. Let us integrate instead Xee( r 5 T = 0) over r first, and then over y since we are 
at T = 0. The result is the static electron-electron contribution to the succeptibility : 

Av (T=0) = J- fan \S{E,T B )-g(EM 
X " ( ' « E2 ' (5.38) 



4n 2 T B ' 

where now the integration over y made sense since the E integral was convergent. Hence 
we conclude 

A Xee (r,T = 0) = ^-5(r). (5.39) 
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Let us now come to the electron-impurity contribution and show that it is equal and 
opposite to the previous contribution. It needs a little trick to be converted to a current 
correlation: Look at the quantity : 

^^[SL P (y') - s! mp (-oc)\. (5.40) 

For finite y, the term < <9 r </>(r, y)Sf mp (— oo) > will be zero because of the infinite time 
separation between the operators, and the average of ( p.40| ) will reproduce what we want. 
On the other hand, the difference of spins can be computed following the transformations 
in section 5.1. We then find that : 

1 f 

- < 3 r 0(r, y)S z imp {y') >= -1- lim / d£ <: d r </>(r, y)dx<t>{X, £) :> . (5.41) 

2 y Sng x^0j_ oo 

In order for the integral over £ to converge, we need to take now the modular transformed 
representation of the current correlation. We finally obtain the expression : 



-i- I l§[F(E,T B )-F(E,0))e- iEr + E y' + c.c. 
8^9 Jq E 



(5.42) 



87TQ 

Then using analytical continuation, we find the contribution to be the same as the previous 



one up to a sign Ax e i(f) = — Ax ee (r). This is consistent with Lowenstein [32 1 . This is 



also consistent with results recently found in perturbation theory at finite T |33[ showing 
that these succeptibilities are zero everywhere outside r = 0. This is inconsistent with the 
results of []33| in that in their case, the extrapolation of the perturbative results towards 
zero temperature show that both Axi e and Ax ee go to zero (not the sum). This last 
behaviour is also expected from a physical argument. We don't have a clear understanding 
of the discrepancy yet, since making contact between the two calculations is not easy. Our 
guess is that there is a subtlety when the T — > limit is taken at the same time as the 
coupling to the impurity and g — > 1, and that working directly at zero temperature is not 
the correct way to proceed. We hope to be able to verify this explicitely by making a finite 
temperature calculation in the same formalism soon and check whether the discrepancy 
comes from the limit g — > 1 or T — > 0. In that sense, the last results, should be taken with 
caution. 

All these results are related to the so called uniform part of the succeptibility. The 2k f 
part, defined in section 2, involves a different set of operators once bosonised: e^l 2 . These 
operators have differnet anomalous dimension in the IR and UV and their treatment using 
massless form factors encounters serious complications which are discussed in the next 
section. 
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6. Form- factors failures 

For operators with naive engineering dimension such as the current, the massless 
scattering approach is thus seen to work quite nicely. Another candidate for which things 
work is the stress energy tensor, that is basically (d x (j)(x)) , describing the density of energy 
at some distance from the boundary (the impurity). Unfortunately, things are very different 
for operators with a dimension that is not constrained by any symmetry. An example is 
cos h(j){x). This operator is of crucial physical interest: for instance knowing its correlators 



would lead to an exact determination of Friedel's oscillations fl34|. Unfortunately key 



difficulties appear here with the form- factors approach. Form- factors of the operator cos hi 



in the bulk massive sine-Gordon model have been determined ||15|| . A first feature is that 
they depend only on rapidity differences, and do not exhibit any overall factor depending 
on the rapidity scale and fixing the (naive) dimension, as for the current. We can then 
take massless limit as described above. What one finds is that, at least for q = t— r- — , the 

' f integer ' 

only form- factors whose limit is not cancelled exponentially by powers of e~^° are those 
which are left and right neutral - for instance form-factors between the ground state and 
any state made only of breathers. A first problem then arises when one tries to reproduce 
properties of the operator cos h<p in the bulk massless theory. For instance one finds that 



the one breather form-factor [35]is a pure, rapidity independent, number, 



< 0| cos —\9 >i= c 

In the massless limit, its contribution to the two point function will therefore be of the 
form 

/oo 
e Mze "d(3 
-oo 

an integral that is IR divergent - recall that for the current it was the term giving the 
naive engineering dimension that made the integrals converge in the IR. One might think 
that this problem has a simple solution: put by hand the anomalous dimension, ie multiply 
in the massless limit each of the form-factors by a factor m 9 / 2 , where m is the bulk mass. 
The above contribution reads then 



/oo 
exp m [cosh (3x + i sinh (3y] d(3 
-oo 



Unfortunately, in the limit m — ► 0, it does not behave any better, and does not reproduce 
the expected \x\~ 9 behaviour. Recall how, for the current, the dimension followed trivially 
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from the scaling properties of each individual form-factor, and was reproduced by every 
term. Here, the 1-breather example suggests that the correct behaviour can be reproduced 
only when the whole series is summed up at finite m and then only m is sent to zero. 
There is thus no hope to reproduce even approximately the bulk-behaviour by massless 
form- factors. A related situation has been encountered in massless flows in ]13| . 



We have tried to study the effect of the boundary by putting in by hand the correct 
massless bulk behaviour and concentrating on the corrections induced by the boundary - 
for instance by studying ratios of correlators with different boundary couplings. Similar 
problems unfortunately occur. A simple way of seing this is to consider the expected form 
of the correlations of this operator in the UV and IR limit. In the IR the field <p obeys 
Dirichlet boundary conditions, and one has 

/ 9 

I -I 7.-7 

< 4>{z,z)(j){z ,z >= — 4(7 In 



so, restricting to arguments on the x axis, 



< cos — (x) cos -{x ) >- 



z + z' 



x + x' 



(6.1) 



x — x' 

In the UV on the other hand it obeys Neumann boundary conditions and thus 

< cos —(x) COS —(X ) >- 
2 2 



_ T '\ 2 9 ' 



\X — X 



(6.2) 



(6.3) 



We see thus that the anomalous dimension of the field cos % is changing, being equal to g/2 
in the IR and to g in the UV. This is of course very different from the case of derivatives 
of 4>. Thus, at the present time, the form-factors approach fails for that operator. 



7. Conclusions. 

In this paper we presented a method to obtain current-current correlations in massless 
theories with interaction at the boundary. The technique, as presented, is quite general 
and has been applied to different problems successfully. Several generalizations along the 
same lines appear possible. The most interesting would be to study similar quantities but 
in the presence of a bias (or voltage) and temperature. In that case the ground state is no 
longer empty but filled with quasi-particles in a way determined by the thermodynamic 
Bethe ansatz. A combination of the TBA and the technique presented here can hopefully 
allow the determination of the current-current correlation in that case. 
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Other generalizations appear more challenging. As explained in the previous section, 
the approach, when applied to operators with a non trivial dimension, naively fails, pre- 
venting us eg to study Friedel oscillations for the moment. The approach also fails when 
one gets too close to the isotropic point g = 1. Although this should not be too catas- 
trophic in practice because results can be extrapolated from the g < 1 regime, this remains 
a big challenge. 
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